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' Abstract 

00 

We studied the possibility of the spin pump in a S=l/2 antifcrromagnetic 
chain. The spin chain is mapped into a fermion system and bosonization is 
utilized to transform the equation of motion to a sine-Gordon equation. The 
D I sine-Gordon equation on a finite chain with different boundary conditions is 

^ • solved. Among numerous solutions, the static soliton is compatible with the 

^ ' original physical system. By varying adiabatically a angle (p in the phase space 

composed of applied electric and magnetic fields, the spin states change between 
^ I the Neel state and dimer state and a quantized spin S* = 1 is transported by 

^ i the bulk state from one end of the system to the other. 

g . PACS: 75.10.Pq, 75.lG.Jm, G3.65.Vf 

Q ' Keywords: spin chain, sine-Gordon equation, quantum spin transport 

(N 
> 

in 

00 

(N 

G\ 
O 
oo 
O 

> 



1 



1. Introduction: An adiabatic quantum pump is a device that generates a 
dc current by a cyclic variation of some system parameters, the variation being 
slow enough so that the system remains close to the ground state through- 
out the pumping cycle. After the pioneering work of Thouless[l] and Niu and 
Thoulesss[2], the quantum adiabatic pumping physics gets more attention. It is 
applied to the systems like open quantum dots [3-5], superconducting quantum 
wires [6, 7], the Luttinger quantum wire [8], the interacting quantum wire [9] and 
of course the spin sytems. 

In recent years, spintronics become an exciting new field new field of re- 
search. Various proposals of generating spin current have been studied. Among 
them, an adiabatic spin pumping process is most interesting. Quantum spin 
pumping physics probably has inspired by the phenomenal work of Thouless[l], 
which is clearly related to the topological explanation of quantum hall effect by 
Thouless et.al.[10]. However Halperin[ll] pointed out before that the quasi-one 
dimensional edge states played an important role in quantum hall effect. Hat- 
sugai[12,13] showed that the edge states indeed have topological meaning and 
thus confirm their importance. 

Shindou[14] has shown that the origin of spin transport is due to the edge 
state of the system. Fu, Kane and Mele[15] and Fu and Kane[16] studied the 
similar problem. Among their contribution, they found that the edge state 
crossing (Kramers degeneracy) is essential for spin pumping. So the possibility 
of spin transport through the bulk states of the system is not reveal from these 
studies and leave this bulk state spin transport as a open problem. 

Here we mention very briefly the basic theme of spin transport in adiabatic 
process. Suppose we consider a spin chain and constructed a parameter space 
with (/ist,A) = R{cos(j),sm(j)) where hgt is the applied magnetic field A the 
dimer states bond strength. Fixing R and varying (j) adiabatically in time, one 
can argue that a line integral of A„(K) = {i/2-K) < n(K)|VK|'^(K) >, where 
n(K) is the Bloch function for the n-th band, and K = (fc, A, /;.,,/), on a closed 
loop yields exactly ±1 due to the singularity at the origin. In other words, 
A is related to a fictitious magnetic field B„(K) = Vk x A„(K). One with 
the Stokes' theorem, can express the line integral in terms of surface integral 
(/ B • dS) where the integration is on two dimensional closed surface enclosing 
the origin. This is exactly the quantization of particle transport proposed by 
Thouless [I]. It is well known to us that one can express spin chain problem 
into a spin less fermion problem with Jordan- Wigner transformations. One can 
use this kind of adiabatic variation of parameters as a tool of quantized spin 
transport. Shindou considered the spin polarization Pgz = X^jLi i^j^' ^''^'^ 
divided it into two parts, bulk state part and edge state part and he concluded 
that edge state part of spin polarization contribute to spin transport. Fu and 
Kane[16] considered a similar system with an additional interaction of spin-orbit 
coupling. They calculated the energy bands of the bulk states and end (edge) 
states and were able to show clearly that whenever there is Kramers degeneracy 
of end states, there is spin transport and it has a symmetry. 



2 



We have already mentioned that in ah previous studies of spin transport, 
the contribution is coming from the edge states. Here, most probably first in 
the literature, wc raise the question, whether the edge states arc indispensable 
in spin transport?. We do the rigorous analytical exercises to complete the 
search of this question. One can see during our analytical derivation that spin 
transport is nothing but the transport of soliton in the system. The plan of 
our paper is the following: In section (2), we present model Hamiltonian and 
continuum field theoretical studies. Section (3) is for analysis of sine-Gordon 
equation on a finite chain for different boundary conditions. Section (4) is for 
the detail analysis for static soliton solution. Section (5) is reserved for results 
and discussions. 

2. Hamiltonian and Continuum field theoretical studies We consider 
a spin 1/2 chain of finite length, described by the Heisenberg Hamiltonian sim- 
ilar to that of Shindou[14]. A controlled dimcrization amplitude and applied 
magnetic field are also present. The total Hamiltonian has three parts: 

H{t) =Ho + ifdim + H,t (1) 

where 

N 

ifo = J^S,-S,+i, (2) 

i=l 

iJdim = ^ + 5r5+ J, (3) 

and 

H,, = h,,{t)Y^{-iys!. (4) 

i 

Hdim is the bond alternation term which can be induced by applying an electric 
field to the spin chain to alter the exchange interaction. It introduces into the 
system the strength of dimerization A{t). Hgt is the coupling of the system to 
a staggered external field hgtit). The time-dependent bond strength A(t) and 
staggered field hst{t) can be varied adiabatically so as to create a parameter 
space for Berry phase. We write A and hgt as A) = R{cos(f),sm(f)), with 
R fixed. Varying (f) adiabatically, we expect spins to be transported. We shall 
also argue that going through one cycle along the loop, there will be quantized 
spin component transported from one end to the other. 

The method of bosonization [17-20] has been used successfully to treat various 
one-dimensional systems, including the spin chains. It is suitable for the system 
we are considering. To this end, we first make the Jordon- Winger transformation 
to represent spins by fermion field fi and . Then, the bosonizations of fi and 

will performed. 

i-i 

/,=exp(i7r^5,+5-)5r, (5) 
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and 



where 



i-i 

/; = 5+expH7r^5+5r), (6) 

fj - R{xj)e''"'''J + L{xj)e-'^'"'i , (7) 

/t - i?t(a;,)e-''=^^^- + L\xj)e'''^='^ , (8) 



and 



R(^x) = ^Ur/ie^t^+(^)+^-(")l/2 (9) 
V27rQ; 

L{x) = -^Tj^e'^-^+(-)+^-i-)]/^ (10) 

are the slowly varying fields, and r]i and r/2 are the Klein factors. Here, kp is 
the Fermi wave vector and and a is the lattice constant. For a half filled system 
we have kp = 7r/2a. In the following derivations, we left out the details because 
they can be found in many textbooks[17-20]. 

Sj = fth - \ (11) 



'3 Jj Ji 2 

= ^|^-(-lF;^sin^,(.,) 

and 



^t^i+l + — fj^fj+i + fj+lfj 

= -a[ATrfP + -^{dj+f] - i-iy— cose+{xj). (12) 

Here 0+ = is bosonization phase and Il{x) = — {1/4^)8x0 -{x) is the conjugate 
momentum of ^+(x). Substituting eqs. (11) and (12) into eqs. (1-4), and drop- 
ping the rapidly varying components such as !)•' cos^+(a;j), we obtained 

j 

H = [ dx{v[-Kr]lP + -^{dj+f] - sin(?+ + ip) + cos2?+} (13) 

where the velocity 

and the quantum parameter 

V = ^ (15) 
were discussed in ref. 17-20. Thus, we have the equation of motion 

0^0+ = v'^dl0+ + cos(?+ ^"f)^^ sin 20+ (16) 
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The term of sin 2(?+ is irrelevant in the sense of renormahzation group analysis, 
so we consider only the part 

d^e+ = dle+ + cos{e+ + ^p). (i?) 

where we have change variables: z = \/2JRx/ua and r = ^/^IjRt/a. It is 
similar to the standard sine-Gordon equation 

dle+ - dle+ + sin 6>+ = (18) 

which has been well-studied. Howevc;r, for our purpose which is to stiidy the 
spin transport, we will solve it on a chain of finite length on which the phase </> 
is no longer a trivial constant but introduces new meaning to the solution. This 
way, one can recognize the motion of spins from one end to the other. 



3. Analysis of sine-Gordon equation on a finite chain We shall analyze 
eq. (18) first. The result can be applied to eq. (17). Equation (18) has 
many kinds of solutions. The traveling-wave solutions, such as arctan[exp(7(z — 
vt))] is not suitable for our purpose because they cannot meet fixed boundary 
conditions. For the finite- length systems, we consider the so called "separable 
solutions" [21-23]. 

(A(^,t) =4arctan(^^^). (19) 
and /(/3-2;) and giP-r) must satisfy the following equations:. 

id.ff = + M/' + (^)] (20) 

and 

i^rgf = (jj2)[- V + (m - 1)5' + «]• (21) 

with the requirements ^j?+Ak\ > and (^— 1)^-|-4kA > 0. A, ^, k, X f3 and fl are 
mutually related constants. We will show how they are determined in a while. 
First, we would like to put forward the observation that f{f3z) and g{^lT) satis- 
fying eqs. (20) and (21) are Jacobi elliptic functions (JEF)[24]. Jacobi elliptic 
functions are defined as the following: 



sn{u) 1, 

''^ . (22) 

V(l-t2)(l- A;2f2) ^ > 



where sn{u) is one of the JEFs and fc is a constant in the range [0,1]. The 

second JEF is cn{u) where sn'^{u) + cn'^{u) = 1. There are more JEFs. They 
can be found in Appendix A. The ones we are going to encounter are sc(tt) = 
sn{u)/cn{u) and d'nP{u) = 1 — k'^s'n?{u). Both / and g are JEFs and their 
constants are denoted by kj and kg. fi, k and A are constants determined by kf 
and kg. The relations are different for different Jacobi elliptic functions. 
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Here wc give an example of f{l3z) = cn{j3{z — 2:0)) and g{^T) = cn(f7T) 
where zq is a constant. With the equation for cn(u) (see Table I in Appendix 
A) 

{duCn{u)f = (1 - w2)(l _ fc2 + )t2y2-)^ (23) 

we found from comparison with eq. (20) that kA^ = kj^, fj, = 2fcj — 1 and 
X/A^ = 1 — k'j. As a result, we get 



^2 ^2 



and, 



where 



^ 



If we choose the fixed boundary condition 9+{z = 0) = 6+{z = L) = with L 
being the length of the system, then we will have 

(3L = AlK{kf) (27) 

with 

K=[\ (28) 

io ^{l-t^){l-kH^) 

being the complete elliptic integral of the first kind, zq = L/Al and I is an 
integer. 

Not all the combinations of JEFs can satisfy the sine-Gordon equation. A 
table of the differential equations for all the JEFs is given in Appendix A. We 
will discuss the solutions of the sine-Gordon equation in eq. (17) on a finite 
system under various boundary conditions. Although finite-length solutions are 
well-known, different boundary conditions and the presence of (j) will impose 
restrictions on the solutions and infuse physical meaning to the wave forms. 



Case 1: Periodic boundary condition 9+{z,t) = 9+{z + L,t) The 
first boundary condition coming to mind is the periodic boundary condition. 
There are many combinations of JEFs that can satisfy the periodic boundary 
condition. Here are two examples. 

la: 

9+{z, t) — — — (fi -|- 4arctan{Acn[/3(z — zq); kf]cn[nT; kg]} (29) 

where /3L = AlK(kf) and 
lb: 

^ TT 

9+{z,t) = — — (fi -\- 'iaxctaxi{Asc{j3{z — ZQ);kf\dn^T;kg\} (30) 
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where /3L = 2lK{kf) and zq is a arbitrary constant. For this boundary con- 
dition, there is no spin transport if one varies the parameter (p adiabatically. 
The reason is quite simple. Increasing only give a constant change to 0+{z) 
everywhere. Thus the fermion field operators on every site from Jordan- Wigner 
transformation acquire a constant phase and the spins remain the same. 

Case 2: Fixed boundary condition 0+{z = 0,t) = d+{z = L,t) = It 
seems that we can have the solutions like 

^ TT 

6+{z,t) = — — if + 4 arctan{ Acn[/3(z — zo); kf]cn[flT; kg]} (31) 

where /3i = 2lK{kf). However, the presence of the adiabatic change term 
7r/2 — ip in front requires that cn[/3(z — za); kj] to be finite. Then the function 
cn[ftT; kg] makes the inverse tangent function varying with time and hence, the 
forms in solution (32) can not satisfy the fixed boundary condition. 

Case 3: Free end boundary condition {d9+{z, T)/dz)]z=o = (30+ (z, T)/dz)]z=L = 
The solution is 

0+{z, t) — — — if -I- 4arctan{A(in[/3(z — zq); kf]sn[QT; kg]} (32) 
where l3L = 2Kf and (izo = Kf. The energy is equal to 16PE{K) where 

7r/2 


is the complete elliptic integral of the second kind. This solution cannot provide 
the system with spin transport for the same reason as that in case 1. 



4. Detailed analysis of the static soliton case It is most interesting to 
study the solution of the static soliton of eq. (17). Let us first consider the 
boundary condition 

^+(^ = 0) = 0, (34) 
e+{z = L) = 27r. 

The phase difference 27r implies that the fermion field or the spins have same 
boundary conditions at both ends and h{ui('(\ a common case for a finite spin 
chain. On the other hand, it is a fixed boundary condition of 0+. Therefore, 
different values of (j) will induce distinct solutions. The soliton has the form 

^ TT 

0+(z,t) — — — LP + 'iaict&n{Asc[l3{z — ZQ);kf]dn^T]kg]} (35) 

where 
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n = AP, (38) 

and 

I3L = K{kf). (39) 

Equations (36-38) can be derived by substituting eq. (35) into cq. (17) and 
eq. (39) comes from the boundary conditions in eq. (34). We seek the static 
solution because it can always satisfy above boundary conditions. In this case, 
we require kg = 0, dn{Q,T, kg = 0) = 1 and A takes a special value Ath- In view 
of eqs. (36) and (37), the static soliton is 

0+{z,t) = - </5 + 4:a.vcta.n{AthSc[(3{z - zo);kf]} (40) 

with 

/3 = 1/(1-A?J (41) 

and 

kf = ^l-^t. (42) 
The boundary condition at ^; = requires that 

tan(^ -^)= AthSc{-(}zo) (43) 

which determines zq. For the boundary condition aX z = L, we have derived the 
following lemma: 

For the solution in eq. (40) with 



0L = m, = ^^>=J^ 



dt 



V(l-i2)[l-(l-A^,)i2]' 



(44) 



the difference of 0+{z = 0) and 0+{z = L) is always equal to 2-k. 

The derivation is given in Appendix B. Eqs. (40) and (44) are the main 
result of this paper. Equation (44) can be generalized as 

(3L = IK (45) 

where I is any nonzero integer. As The larger the magnitude of I, the higher 
the energy. 

In Fig. 1 Atfi evaluated with eq. (44) is plotted. It shows that At^ decreases 
rapidly with increasing L. The magnitude of Ath is closely related to the wave 
form of 9+. A small At^ results in a steep change in 9+., or a sharp domain 
wall. As it will be shown later, it related to the quantum spin transport. 

The energy of the static soliton can be calculated with 
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It can be shown that 

1 F)9 1 r)G 

-(^)' +sm(«+ +^) - -{-^)%„_^ - 8m(9+ +0)|„„ (47) 
Hence, eq. (45) becomes 

0/12 ^ 

£ = dz[l+ (Y^^ - 2 sin(^+ + </,)] (48) 

where the term of the time derivative is dropped for we are considering the 
static case. We can change the variable of integration and get 

£ = V2i de+[l + j^, - sin(^+ + 0)]V2 _ L[l + ] (49) 

where 0+^1 = 7r/2— (/?+4 arctan(yl(;isc(— /3zo)) and 0+^2 = 7r/2— (^+4 arctan(ylt^sc(/3i— 
/32:o)). It has been shown in Appendix B that ^+,2— ^+,1 = 27r. Therefore the 
total energy £ is independent of zq and (/) because the integration is over an 
entire period. £ depends on only one parameter, /3, for static soliton because 
A = Ath and Ath is also determined by (3. The spectrum is plotted in Fig. 2 
with eq. (45) . It is very similar to that of standing wave with (i being the wave 
vector. 

In the limit L ^ 00, Ath becomes vanishingly small as it can be seen from 
eqs. (44) and (41). We thus have P, kf ^ 1 and K{kf) ~ \n{A/Ajf^). With eq. 
(41) we get 

^th^2exp(-L/2). (50) 

The magnitude of Ath is small even for a modest length of L. For example,when 
L = 24, Ath ~ 2/e^^ ~ 1.23 X 10~^. Since Ath can be viewed as the amplitude of 

the nonlinear wave, the wave form on a long spin chain becomes flat everywhere 
except for narrow regions cn{z — zo) ~ 0. Therefore, one can expect a sudden 
change of 9+ or a sharp domain wall. 

5. Results and discussion In this section, we present the results of our 
calculation above. First of all, we plotted zq versus <p in Fig. 3. zq can be 
viewed as a reference point of the solution of the sine-Gordon equation. Its 
movement is a clear indication that the solitons is set in motion by (j>. Its 
motion is not smooth as one can easily see that there is an abrupt change near 
= 7r/2, a manifestation of the nonlinearity of the solution. It decreases by a 
distance L when increases by 27r. To see more clearly how the soliton moves, 
we plotted in Figs. 4 9+ versus lattice sites for different values of at L = 24. 
We can go back to the original spin system to see how spins are transported by 
utilizing eqs. (11) and (12). Hence, Sj and Sj^^iSj + S'^^j.^S'^ are also plotted. 
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In view of cq. (11), the "domain wall" or the region where there is a jump of 
0+ is the place where (S^) is large. Hence, in Figs. 4 the jump of 0+ and the peak 
of {S^') move together with varying (f> and spins moves from right to left. These 
figures also show that the static soliton solution really is an Neel state in the 
spin chain except in the neighborhood of </> ~ n/2. In this range, the Neel state 
becomes unstable due to the dimer coupling A{t) we added. This is manifest in 
Figs. 4(b-d) where S'^-^S'" + Sj^iS^' which is proportional to the dimer state 
amplitude is large. Recall that 4> is defined in {hst,^) = -R(cos(/>, sint/)). The 
dimer strength A is the the largest when ^ ~ 7r/2. It is when the Neel state 
becomes unstable and the dimer state amplitude becomes significant that the 
transport of spin becomes possible. Not coincidentally, one can find in Fig. 3 
that zo changes abruptly in this range. 

In Fig. 5, Sj and S^^^-^^Sj^ + S'^^j^S'^ versus lattice sites for = 7r/2 at a 
shorter length of L {L = 14) is plotted for comparison. For smaller L, the curve 
of 6+ is smoother or the domain wall is not as sharp. On the other hand, the 
edge (end) effect is more important. The directions of spins are less ordered for 
a shorter spin chain because the edge effect penetrates deeper into the "bulk". 

We will elaborate more on how the spins are transported. This can be done 
with eq. (11). The spin polarization is 

L 

Ps' = jJ zS'{z)dz. (51) 



By integration by parts, we found 

L L z 



Pg. = J S%z')dz' -j^jj S\z')dz'dz 





(52) 



To find the variation of Pg^ due to (f) we note that the first term remain constant 
as (j) varies. This can be seen by substituting eq. (11) into the integration. The 
contribution of the oscillatory term vanishes as L ^ oo and the term of the 
derivative gives unity due to our boundary condition, no matter what the value 
of (j) is- Thus, denoting the variation of Ps' due to the adiabatic change of (j) by 
6Ps' , we have 

L z L z 

SPs^ = -^{[ J J S^{z')dz'dz]\^=^, -[j j S\z')dz'dz]\4.^^,} 



L L 

^ -^{[ / 0+{z)dz]U. - [ j e+{z)dz]U,} (53) 



where the second step can be reached by using the approximation 

~ {de+/dz)/2-K (54) 
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in integration for largo L. In view of Figs. 4(b-d) where (f) increases exceeding 
7r/2, we found that 0+ increases in the entire length of the system by approxi- 
mately 27r and hence, dPs^ — —1 around = 7r/2 and a spin 1 is moved from 
right to left around ip = 7r/2. In Figs. 4(e-h), where 9+ is almost constant away 
from ends, we did not see any spin movement in the bulk but rather, there are 
changes of spins at both ends. 

To see the quantum spin transport (a spin of unity being transported) more 
clearly, we consider the limit L ^ oo which can be simulate very closely by the 
case L = 24. The following equation will be very useful for our purpose 

06+ ^ AAthPdnjPjz - zp)) 

dz cn^{l3{z-zp)) + Al^sn^{l3{z-zp)y ^ ' 

since it is the dominant contribution to . From eqs. (41), (50), (A-11), (A-17) 
and (A-18), we find that 

de+ ^ 8e-^/^ cosh(z - zq) 
dz ~ l + 4e-'^sinh2(^-^o) 

The peak of d9+ / dz or is at 

z-zqc^L/2. (57) 

The larger L, the narrower the peak. On the other hand, zq is determined by 
the boundary condition. When = ■7r/2 — 5 where (5 is a small and positive 
number, we have 

t&n{--) = AthSc{-l3zo) (58) 

and we find that zp ~ L/2 as long as 5 remains finite (see eq. (50)). The 
resulting due to eqs. (11) and (55) has a peak at the right end and vanishes 
everywhere else. When cj) = 7r/2, eq. (57) gives = and the peak of S'^ moves 
to the center of the spin chain. If ^ = it/2 + 6, then zq ~ —L/2, and the peak 
moves to the left end. 

As for the quantity of spin transported, we can analyze the variation of 9+. 
In view of eq. (40), as cj) = ■k/2 — S, 9+ almost vanishes for the entire chain 
except for the right end. As (j) increases to 7r/2 + 5, ^+ ~ 27r for the entire 
chain except for the left end where it drops to zero sharply. Hence, according 
to eq. (53), during the interval = 7r/2 — S to (j) = tt/2 + 6, a spin of unity 
is transported from the right end to the left end. As </> increases onward from 
tt/2 + 6, there is no spin transport in the bulk. Nevertheless, the plateau of 9+ 
is lowered (see Figs. 4(e-h)) and as a result, at the left end decreases and a 
peak of at the right end start to grow. This kind of change continues until 4> 
reaches 57r/2 — 5. At this stage the state of soliton returns to that of = 7r/2 — (5. 
We conclude the analysis of our result by the following summary: There is a 
swift spin transport in the bulk during the short interval between (j) = it /2 — 5 
and (j) = it/ 2 + 5 where 5 can be made arbitrarily small if L ^ oo. The net spin 
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transported is unity. Beyond this interval, the spins at both ends vary with (f) 
but there is no spin change in the bulk. 

The question will inevitably be raised: Can spin be transported? We have 
seen that the sohton returns to the starting state if cf) increases by 27r. Thus 
there is no net spin transported in a cycle. However, in a realistic system, two 
ends of the spin chain must be connected to leads. The leads ought to serve as 
a spin source and a spin drain. Thus it is reasonable to envisage the following 
picture: At </> = 7r/2 + (5 the left end can dump spin into a spin drain and 
the right end can extract spin from the source. When <f) = 5^/2 — 5 which is 
equivalent to = 7r/2 — (5, the dumping of spin at the left end is complete and 
the peak of spin at the right end has grown into saturation. Then an unity of 
spin is transported from the right end to the left end when (f) increases from 
(/) = 7r/2 — (5 to (/) = 7r/2 + 5. This is in all intent and purpose, same as a physical 
system of spin transport. On this point, our system is same as the Z2 spin 
pump proposed by Fu and Kane[16]. However, there is an important difference. 
Our spin transport is through a bulk state. This is completely different from 
Shindou's[14] and Fu and Kane's[15,16] pictures in which the level crossing of 
the end (edge) states is essential. Consequently, our spin transport is quantized 
because it is through a bulk state. Connecting to spin reservoir cannot destroy 
the quantization as it will do the transport due to end states. 

This work is supported in part by NSC of Taiwan, ROC under the contract 
number NSC 95-21 12-M-002-048-MY3. 

Appendix A: In this appendix we listed some properties of the Jacobian 
elliptic functions. See eqs. (24-25) for the definition. 



sn^('u) + cn^(u) 
dn^{u) + k'^sn'^{u) 
sn{u + K) 
cn{u + K) 



1, 

cn{u)/dn{u) 
—sn{u)/dn(u) 



1 



(A-1) 
(A-2) 

(A-3) 
(A-4) 



where fc'^ = 1 — fc^. The derivatives of Jacobian elliptic functions are 



9„sn(u) = cn{u)dn{u), 
ducn{u) = —sn{u)dn{u), 
dudn{u) = —k^sn{u)cn{u). 



(A-5) 
(A-6) 
(A-7) 
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Using above equations we found the differential equations to be satisfied by the 
Jacobian elhptic functions and listed them in Table 1 where 



ns(u] 


= l/sn{u), 


(A-8) 


nc{u) 


= l/cn{u), 


(A-9) 


nd{u) 


— l/dn{u), 


(A-10) 


sc{u) 


= sn{u)/cn{u), 


(A-11) 


sd{u) 


= sn{u)/dn{u), 


(A-12) 


cd{u) 


= cn{u)/dn{u), 


(A-13) 


cs{u) 


= l/sc{u), 


(A-14) 


ds{u) 


— l/sd{u), 


(A-15) 


dc{u) 


= l/cd{u). 


(A-16) 



Having checked those equations in Table I, one can see that there are many 
combinations of JEF that can satisfy eqs. (20) and (21) where the discriminant 

is n"^ + 4fi;A for kf and {fi + 1)^ + 4kA for kg. 

For L — > oo, we find from eq. (44) that k ^ 1 and 

sn{u) ~ tanh(M), (A-17) 

cn(u) ~ rfn(u) ~ sec /i(w), (A- 18) 

Appendix B In this Appendix, we prove the following lemma. 

Lemma 1 // f?+(z) = 7r/2 + + 4arctan{At/iSc[/3(L — ^o); A/]} andA = Ath, 
then 0+{z = L) — 6+{z = 0) = 21-k where I is a natural number defined by 

withp = l/{l-Al). 
Proof: 

The boundary condition requires that at ^ = 

^ = — arctan^(/isn(— /32:o) (B-2) 

and at z = L 

TITT (b 

— — — = aictan AthSc{P{L — zq)) = aictan AthSc{lK — (3zo) (B-3) 

Note that the period of sc{u) is 2K. Hence, when z increases from to L, 1/2 
periods will pass. Since the period of arctan function is tt, the term 1^/2 has to 
be added on the left hand side of eq. (C-3). To see it more explicitly, define 

a = aTctan{Ath,sc{(3{L — zq))) = aictan{AthSc{lK — (3zq)), (B-4) 
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and 

^ = arctan(At/iSc(-/3zo)), (B-5) 

then we have 

tana- tan ^ sc{lK - ^zq) - sc(-(3zo) 

tan(a — <;) = = At^ — - — - — - — -. (B-6) 

1 + tan a tan 1^ 1 + Afj^scilK — l3zQ)sc[—pzo) 

With eqs. (A-3) and (A-4), we have 

sn{nK - Pzg) cnjnK - K - Pzq)) 
sciPinK - zo)) = = -k'^sninK - K - Pzo))- ^^"^^ 

where kj + kf = 1. Since kj = 1 - A^^, we find that fc^ = Af^. Eq. (C-6) 
bcconiGS 

csaif-ir-/3zo)/fc}-=sc(-/3zo) 

tan(a - ,) - "-^^^ I - sc{-Pzo)csilK ^ K - Pz,) (^"'^ 

where cs(u) = l/sc{u). So the denominator of tan(ck — <;) vanishes as n = 1 
and the numerator is finite. Thus tan(a — <;) = ±oo and a — <; = ±7r/2. The 
sign is determined by the fact that as z increases, 9^ also increases, thus the 
positive sign should be chosen and 0+{z = L) — 9+{z = 0) = 2tt. If Z = 2, we 
can use eqs. (A-3) and (A-4) again or simply use the fact that the period of 
sc{u) is 2K to find out that the numerator vanishes. Thus, tan(Q! — <;) = tt and 
6+{z = L) — d^{z = Q) = Att. Hence, we conclude that 0^{z — L) — 9^{z = 0) = 
4(a - <?) = 21tt. 
End of Proof. 



References: [1]. D. J. Thouless, Phy. Rev. B27, 6083 (1983). 
[2]. Q. Niu and D. J. Thouless, J. Phys. A, 17, 2453 (1984). 
[3]. P. W. Brouwer, Phys. Rev. B58, 10135 (1998). 

[4]. T. A. Shutenko, I. L. Aleiner and B. L. Altshuler, Phys. Rev. B61, 
10366 (2000). 

[5]. Y. Levinson, Entin-O.Wohlman and P. Wolfe, Physica A302, 335 (2001); 
Wohlman-O.Entin and A. Aharony, Phys. Rev. B66, 35329 (2002). 
[6]. M. Blaauboer, Phys. Rev. B65, 235318 (2002). 

[7]. J. Wang and B. Wang, Phys. Rev. B65, 153311 (2002); B. Wang and 
J. Wang, Phys. Rev. B66, 201305 (2002). 

[8]. P. Sharma and C. Chamon, Phys. Rev. Lett. 87, 96401 (2001) and 
[cond--ma t/0209291, 

[9]. R. Citro, N. Anderi and Q. Niu, cond-mat/0306181' 

[10]. D. J. Thouless, M. Kohmoto, M. P. Nightingale and M. den Nijs, Phys. 
Rev. Lett. 49, 405 (1982). 

[11]. B. I. Halperin, Phys. Rev. B25, 2185 (1982). 

[12]. Y. Hatsugai, Phys. Rev. Lett. 71, 3697 (1993) 

[13]. Y. Hatsugai, Phys. Rev. B48, 1185 (1993). 

[14]. R. Shindou, J. Phys. Soc.Jpn. 74, 1214 (2005). 



14 



[15]. C. L. Kane and E. J. Mclc , Phys. Rev. Lett. 95, 146802 (2005). 

[16]. L. Fu and C. L. Kane, Phys. Rev. B74, 195312 (2006). 

[17]. E. Pradkin, 'Field Theories of Condensed Matter Systems', Addison- 
Wesley, Redwood City, CA (1991). 

[18]. N. Nagaosa, 'Quantum Field Theory in Stongly Correlated Electronic 
Systems', Springer- Verlag, Berlin (1999). 

[19]. A. O. Gogolin, A. A. Nerscsyan and A. M. Tsvclik, 'Bosonization and 
Strongly Correlated Systems', Cambridge University Press (1998). 

[20]. T. Giamarchi, 'Quantum Physics in One Dimension', Oxford University 
Press (2004). 

[21]. G. Costabile et. al., App. Phys. Lett. 32, 587 (1978). 
[22]. R. M. DeLeonardis et. al., J. App. Phys. 51, 1211(1982). 

[23]. R. M. DeLeonardis ct. al., J. App. phys. 53, 699 (1982). 
[24]. Derek. F. Lawden, 'Elliptic Function and Application', Springer- Verlag 
NY (1989). 



15 



Fig captions: 

Fig. 1 Threshold amphtude (Ath) versus system lengrh L for the static 
soUton solution. 

Fig. 2 Energy spectrum of the static soliton versus n = j3/L with the 
boundary condition 9+{z = Qi) =0 and 9+{z = L) = 2mT when L = 24. 

Fig. 3 zo versus adiabatical paramter (p with L = 24 for static soliton and 
boimdary conditions 6'+(z = 0) = and 0+{z = L) = 2n. 

Fig. 4 Sz and dimer state amphtude S^S~ + S~ versus lattice sites 
with L = 24 and different values of <p (a) (p = 0, (b) </) = 7r/4, (c) (j) = '^/2, (d) 
(f) = 37r/4, (e) = tt, (f) (/) = 57r/4, (g) = 37r/2, (h) = 37r/2. 

Fig. 5 5^ and dimer state amplitude S~^S~ + S~S~^ versus lattice sites 
with L= 14 and (j) = n/2. 
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Fig. 4(a) 
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Fig. 4(b) 
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Fig. 4(c) 
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Fig. 4(e) 
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Fig. 4(f) 
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Fig. 4(g) 



34 



1.0 



0.8 



0.6 



0.4 



0.2 



cp=37i/2 



? ? T ? ? 1 ? T ? 



I' \\ K [j " 

. {],-o -O- O Hf| -o - p- O <^ O Hf[ -o 



0.0 <> 



-0.2 



-0.4 



-0.6 



-0.8 



-1.0 



v.s. Site 

— o— e /2n V.S. Site 

+ 

—°— S*S"+S"S* v.s. Site 



12 

Site 



16 20 



24 



35 



Fig. 4(h) 
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Fig. 5 
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Tabic 1 

Differential equations satisfied by Jaeobian elliptic fiuictions. See Appendix 



A for tin" deliuatioiis of Jacol)iau ellii)tic functions. 



JEF 


its equation 


JEF 


its eqiiation 


y = sn{u) 


{duyr = {I - y'){l - k-'y') 


y = cn{u) 


{duyY = (1 - y'){l + k'y') 


y = dn{u) 


{duvY = {y' - - fc' - y') 


y = ns{u) 


(W = (y'-i)(y'-fc^) 


y ~ nc{u) 


{duyY = {y'-im-li')y' + k'] 


y = nd{u) 


{duyY = {l^y'm-k:^)y'-l] 


y = sc{u) 


(9„y)^ = (y^ + l)(l + fc'V) 


y = cd{u) 


{duyY = iy' - m'y' - 1) 


y = sd{u) 


(5„2/)^ = (l-fc'V)(l + fcV) 


y = cs{u) 


iduyY^i^ + y')ik''' + y') 


y = dc{u) 


{duyf = {k'-y'){l-y') 


y = ds{u) 


{duyY = {y'-k'-'){y' + k-') 
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